Detecting quantum non-Gaussianity via the Wigner function 
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We introduce a family of criteria to detect quantum non-Gaussian states of a harmonic oscillator, that is, 
quantum states that can not be expressed as a convex mixture of Gaussian states. In particular we prove that, for 
convex mixtures of Gaussian states, the value of the Wigner function at the origin of phase space is bounded from 
below by a non-zero positive quantity, which is a function only of the average number of excitations (photons) of 
the state. As a consequence, if this bound is violated then the quantum state must be quantum non-Gaussian. We 
show that this criterion can be further generalized by considering additional Gaussian operations on the state 
under examination. We then apply these criteria to various non-Gaussian states evolving in a noisy Gaussian 
channel, proving that the bounds are violated for high values of losses, and thus also for states characterized by a 
positive Wigner function. 

PACS numbers: 03.65.Ta, 42.50.Dv 



I. INTRODUCTION 

Several criteria to detect non-classicality of quantum states 
of a harmonic oscillator have been introduced, mostly based on 
phase-space distributions [fll-nTI, ordered moments lfT2l - IT4l . 
or on information-theoretic arguments [15-21|. At the same 
time, an ongoing research line addresses the characterization 
of quantum states according to their Gaussian or non-Gaussian 
character ll22lJ3T1l , and a question arises on whether those two 
different hierarchies are somehow linked each other. 

As a matter of fact, if we restrict our attention to pure states, 
Hudson's theorem ll32l[33l establishes that the border between 
Gaussian and non-Gaussian states coincides exactly with the 
one between states with positive and negative Wigner functions. 
However, if we move to mixed states, the situation gets more 
involved. Attempts to extend Hudson's theorem have been 
made, by looking at upper bounds on non-Gaussianity mea- 
sures for mixed states having positive Wigner function lf34l . 
In this framework, by focusing on states with positive Wigner 
function, one can define an additional border between states 
in the Gaussian convex hull and those in the complementary 
set of quantum non-Gaussian states, that is, states that can not 
be expressed as mixtures of Gaussian states. The situation is 
summarized in Fig. [T[ the definition of the Gaussian convex 
hull generalizes the notion of Glauber's non-classicality 135] , 
with coherent states replaced by generic pure Gaussian states, 
i.e. squeezed coherent states. 

Quantum non-Gaussian states with positive Wigner function 
are not useful for quantum computation If36ll37l , and are not 
necessary for entanglement distillation, e.g. the non-Gaussian 
entangled resources used in |38| are mixtures of Gaussian 
states. On the other hand, they are of fundamental interest 
for quantum information processing and quantum optics. In 
particular, since no negativity of the Wigner function can be 
detected for optical losses higher than 50% |39| (or equiva- 
lently, for detector efficiencies below 50%) criteria able to 
detect quantum non-Gaussianity are needed in order to certify 
that a highly non linear process (such as Fock state generation, 
Kerr interaction, photon addition/subtraction operations or con- 
ditional photon number detections) has been implemented in 
a noisy environment, even if no negativity can be observed in 
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FIG. 1. Venn diagram description for continuous-variable quantum 
states with positive Wigner function. The quantum states can be 
divided in two sets: quantum non-Gaussian states and states belonging 
to the Gaussian convex hull. The latter trivally includes (Glauber) 
classical states and Gaussian states. 



the Wigner function. 

Different measures of non-Gaussianity for quantum states 
have been proposed ll22] 424l. but these cannot discriminate be- 
tween quantum non-Gaussian states and mixtures of Gaussian 
states. An experimentally friendly criterion for quantum non- 
Gaussianity, based on photon number probabilities, has been 
introduced [26], and then employed in different experimen- 
tal settings to prove the generation of quantum non-Gaussian 
states, such as heralded single-photon states [27], squeezed 
single-photon states |28 1 and Fock states from a semiconductor 
quantum dot |29|. 

In this paper we introduce a family of criteria which are able 
to detect quantum non-Gaussianity for single-mode quantum 
states of a harmonic oscillator based on the Wigner function. 
As we already pointed out, according to Hudson's theorem, the 
only pure states having a positive Wigner function are Gaussian 
states. One can then wonder if any bound exists on the values 
that the Wigner function of convex mixtures of Gaussian states 



can take. By following this intuition we present several bounds 
on the values of the Wigner function for convex mixtures of 
Gaussian states, consequently defining a class of sufficient 
criteria for quantum non-Gaussianity. 

In the next section we will introduce some notation and 
the preliminary notions needed for the rest of the paper. In 



Sec. Ill we will prove and discuss our Wigner function based 



criteria for quantum non-Gaussianity and in Sec. IV we will 
prove their effectiveness by considering different families of 
non-Gaussian states evolving in a lossy (Gaussian) channel. 
We will conclude the paper in Sec. [Vlwith some remarks. 



II. PRELIMINARY NOTIONS 

Throughout the paper we will use the quantum optical termi- 
nology, where excitations of a quantum harmonic oscillator are 
called photons. All the results can be naturally applied to any 
bosonic continuous-variable (CV) system. We will consider 
a single mode described by a mode operator a, satisfying the 
commutation relation [a, a^ = 1. A quantum state g is fully 
described by its characteristic function 1 39 1 



x[gM^^40D{-f)], 



(1) 



where D{'y) — exp{7a — 7*a' } represents the displacement 
operator In addition, the quantum state g can be fully described 
by the Fourier transform of the characteristic function, /. e. the 
Wigner function [.39i 



W[g]{a) 



<f^ 
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(2) 



A state is defined to be Gaussian if and only if its Wigner 
function (or equivalently its characteristic function) is Gaussian. 
All single-mode Gaussian states can be expressed as 

g^D{a)S{OvpS\OD\a), 

where S{^) = exp {^^(a^')^ — ^^*a^}, is the squeezing 



operator, and vp 



-ISa^a 



/Tr[e- 



-fla^a 



] is a thermal state 



(a, ^ £ (D and /3 > 0). Pure Gaussian states can be written 
as If/'c) = D{a)S{£,)\0), and, according to Hudson's theorem 
132, 33], they are the only pure states having a positive Wigner 
function. Together with that of a Gaussian state, one can define 
the concept of Gaussian map: a quantum (completely positive) 
map is defined Gaussian iff it transforms Gaussian states into 
Gaussian states. All unitary Gaussian maps can be expressed 
as Uq — exp{~iHii]\t}, and they correspond to Hamiltonian 
operators i/bn at most bilinear in the mode operators. Similarly, 
a generic Gaussian map can be decomposed as a Gaussian uni- 
tary acting on the system plus an ancilla (the latter prepared in 
a Gaussian state), followed by partial tracing over the ancillary 
mode |40|. 

Another complete description of a CV quantum state g may 
be given in terms of the so-called P-function P[g]{a) 
defined implicitly via the formula 



g^ / d^aP[g]{a)\a){a\ 



(3) 



where \a) = D{a)\0) represents a coherent state. According 
to Glauber a state is non-classical iff its P-function is not a 
proper probability distribution, e.g. the P-function is more 
singular than a Dirac-delta function. Note that the negativity 
of the Wigner function is a more restrictive definition of non- 
classicality: there exists non-classical states having a positive 
Wigner function (e.g. squeezed states), while all the states hav- 
ing a non-positive Wigner function are non-classical according 
to Glauber 

In a similar spirit as in Glauber's approach to non- 
classicality, in this paper we study the concept of quantum 
non-Gaussian states. These are defined as follows. The Gaus- 
sian convex hull is the set of states 



g £T-L\ g 



dAp(A)|V'G(A))(VG(A)| 



(4) 



where p(A) is a proper probability distribution and \iPq{\)) are 
pure Gaussian states, i.e. squeezed coherent states identified 
by the set of parameters A = {a, ^}. Since Gaussian states do 
not form a convex set, the set in Eq. (J4]) includes states which 
are not Gaussian. Moreover any mixed Gaussian state can be 
written as a weighted sum of pure Gaussian states, and hence 
the set above also includes convex mixtures of mixed Gaussian 
states. 

The definition of quantum non-Gaussianity naturally fol- 
lows: 

Definition. A quantum state g is quantum non-Gaussian iff it 
is not possible to express it as a convex mixture of Gaussian 
states, that is iff g ^ Q. 

As illustrated in Fig. [T] the border here defined dividing 
quantum non-Gaussian states and mixtures of Gaussian states 
falls in between the border dividing classical and non-classical 
states, and the one which divides states with positive and non- 
positive Wigner functions. The importance of such a further 
distinction is evident if we note that all states in Q can be 
prepared through a combination of Gaussian operations and 
classical randomization. On the contrary, \f g <^ G then a 
highly non-linear process (due to a non-Gaussian operation or 
measurement) had necessarily taken place in the generation the 
quantum state g. While the negativity of the Wigner function is 
always sufficient to certify it, more elaborated criteria, as those 
elaborated in this paper, are needed in order to detect such a 
characteristic when quantum states exhibit a positive Wigner 
function. 



III. CRITERIA TO DETECT QUANTUM 
NON-GAUSSIANITY 



In order to find criteria for the detection of quantum non- 
Gaussian states, we follow the intuition given by Hudson's 
theorem for pure Gaussian states. We will focus on lower 
bounds on the values taken by the Wigner function of states 
which belong to the Gaussian convex hull Q. In this section we 
present our main findings as one lemma leading to two final 
propositions and two additional corollaries. The 'quantum 
non-Gaussianity criteria', derived directly from these results, 
are presented at the end of the section. 



Lemma 1 (Lower bound on the Wigner function at the origin 
of phase space for a pure Gaussian state). For any given pure 
single-mode Gaussian state |i/'g). the value of the Wigner func- 
tion at the origin of the phase space is bounded from below 
as 



W[\^JG){^JG\m > -exp{-27i(l + n)}, 



(5) 



where n — {ipQ\a^a\ipQ). 



Proof A generic pure single-mode Gaussian state can be al- 
ways written as JV'g) = -D(q;)S'(^)|0), where a = |a|e'* and 
^ = re"^ (r > 0) are two complex numbers. We can thus write 
the Wigner function evaluated in zero as 



W Ug) ii^gim ^ - exp {-2\a\^ [cosh2r 

TT 

+ cos(26l + (/')sinh2r]}. 



(6) 



Our goal is to minimize the value of the Wigner function or 
equivalently to maximize the function 

f{a, = 2|ap (cosh 2r + cos (29 + 0) sinh 2r) . (7) 

A first maximization is obtained by considering 

29 + (f) = 2kiT with fc e N , (8) 

which yields 

/(a, < 2|ape2'' = 2nd (2n, + 1 + 2^^,(1 + n,)) . 

(9) 

In the last equation we introduced the displacement and squeez- 
ing photon numbers, Ud ~ |ap and Ug = sinh r, and we used 
the formula arcsinh(a;) = log(x + Vl + x^)- Note that these 
two parameters obey 

n = {il}Q\a} a\ipQ) = n^ + rig, 

where n is the average photon number of the state \iPq). We 
can thus express the right hand side (rhs) of Eq. (|9]) in terms of 
n and Ug, obtaining 

/(a, < 2(n - Ug) {2ng + 1 + 2^ng{l + n,)) . (10) 

For a given average photon number n, the above function is 
maximized with regard to the parameter Ug by choosing 



l + 2n' 



and obtaining 



This leads to 



/(a,0<2n(l + n). 



(11) 



(12) 



W[\4>Q){'4'G\m> -exp{-2n(l + n)}. (13) 

TT 



By looking at the proof, we remark that the bound obtained is 
tight: given a fixed energy n, by choosing the phases according 
to condition ([8| and the squeezing energy according to \\\\ , 
it is always possible to find a family of pure Gaussian states 
saturating the inequality. In particular the maximization ob- 
tained via condition ([8]l simply corresponds, at fixed Ud and 
Ug, to displace the state along the direction of the squeezed 



quadrature. The condition (111 shows that for small values of 



n, the minimum of the Wigner function is obtained by using 
the energy in displacement, while for larger values of n, the 
optimal squeezing fraction n^ tends to an asymptotic value 
n'^r^ = n/2. 

Let us now generalize the bound obtained to a generic convex 
mixture of Gaussian states. 

Proposition 1 (Lower bound on the Wigner function at the 
origin for a convex mixture of Gaussian states). For any single- 
mode quantum state g which belongs to the Gaussian convex 
hull Q, the value of the Wigner function at the origin is bounded 
by 



W[q\{Q) > -exp{-2fi(l 

TT 



i)} 



(14) 



where n = Tr[ga^a]. 



Proof The multi-index A, which labels every Gaussian state 
in the convex mixture |-0g(-^)) = D{a)S{C)\^)^ contains the 
information about the squeezing ^ and displacement a. We 
can then equivalently consider as variables A = {n, n<;, 6, 0}. 
By exploiting the linearity property of the Wigner function we 
obtain 



W[Qm = j dAp(A)W^[|V'G(A))(V'G(A)|](0) 

>- I dX p{\) exp{-2n(l + n)} , (15) 

where inequality (jSl) has been used. By defining 

/■n /'27r /'27r 

p{n) = dug d(P d0p{X), (16) 

Jo Jo Jo 

which is a valid probability distribution with respect to the 
variable n, Eq. ([TSj becomes 

2 f°° 
W[g]{0)>- dnp(n) exp{-2n(l + n)}. (17) 

T^ Jo 



Studying the second derivative of 



Bmm(") = - exp{-2n(l + n)} 



(18) 



we conclude that the function is convex in the whole physical 
region (i.e. n > 0). As a consequence, 

/ dnp{n) -Bmin(n) > B^m ( / dnp{n) n\ = B^\n{n) 

(19) 



where n = J„°° dnp{n) n — TT[ga^a]. From the last inequaUty 
we obtain straightforwardly the thesis 



W[g]{0) > - exp{-2n(l + n)} 

TT 



(20) 



The following proposition generaUzes the bound obtained 

above. 

Proposition 2. For any single-mode quantum state g (z G, 
and for any given Gaussian map £q (or alternatively a convex 
mixture thereof), the following inequality holds 



W[£z{g)m > - exp{-2n£(l + %)} , (21) 

vr 

where n^ = T!T[£Q{g)a) a]. 

Proof Given a quantum state g which can be expressed as a 
mixture of Gaussian state, and a Gaussian map £q (or a convex 
mixture thereof), the output state 



g' ^Ecig) 



(22) 



can still be expressed as a mixture of Gaussian states. As 
a consequence we can apply to the state g' the result in 
Proposition [T] obtaining the thesis. ■ 



Proof The proof follows from Proposition l2] by considering 
the Gaussian map Ecig) — S{£,)gS^{£,). Moreover, since the 
value of the Wigner function at the origin is invariant under 
any squeezing operation, i.e. 



W[SiOQS\mO)^W[gm 



(25) 



one can maximize the rhs of inequality ( [21) with regard to the 
squeezing parameter ^. ■ 



The violation of any of the inequalities presented in the last two 
propositions and two corollaries provides a sufficient condition 
to conclude that a state is quantum non-Gaussian. We formalize 
this by re-expressing the previous results in the form of two 
criteria for the detection of quantum non-Gaussianity. 

Criterion 1. Let us consider a quantum state g and define the 
quantity 



Alb] =l^[(?](0)--exp{-2n(fi + l)} 



(26) 



Then, 



Ai[e] <o => gig, 

that is, g is quantum non-Gaussian. 

Criterion 2. Let us consider a quantum state g, a Gaussian 
map £q (or a convex mixture thereof), and define the quantity 



Proposition [2] leads to two corollaries that will be used in the 
rest of the paper. 

Corollary 1. For any single-mode quantum state g € Q, the 
following inequality holds 

W[g\{p) > - exp{-2fi^(l + np)} , V/3 e € , (23) 

TT 

where Uf, = Tr[gD{l3)a'faD'f {(3)]. 

Proof The proof is straightforward from Proposition l2] with 

the Gaussian map Ecig) — D{—/3)gD^{~f3). We also use the 
property of the Wigner function 

W[gm = W[D\l3)gD{/3)]{0), 

and D^l3) = D{^ 13). ■ 



Corollary 2. For any single-mode quantum state g belonging 
to the Gaussian convex hull Q, the following inequality holds 



W[g]{Qi) > max ( - exp{-2n^(l + n^)} ) , (24) 



A2[g,£G] = W[£c;{g)m - -exp{-2%(7l£ + 1)} . (27) 

TT 



Then, 



3£:g s.t. A2[g,£c] <0^ g^g. 



Typically, Criterion [Tic an be useful to detect quantum non- 
Gaussianity of phase-invariant states having the minimum of 
the Wigner function at the origin of phase space. On the 
other hand. Criterion [2[ is of broader applicability. To give 
two paradigmatic examples, the latter criterion can be useful 
if: (i) the minimum of the Wigner function is far from the 



where n^ — Ti[gS^^)a^ aS {£,)]. 



origin, so that one may be able to violate inequality ( 23 i by 
considering displacement operations; (ii) the state is not phase- 
invariant and presents some squeezing, and thus one may be 
able to violate inequality ( [24| by using single-mode squeezing 
operations. 



IV. VIOLATION OF THE CRITERIA FOR 

NON-GAUSSIAN STATES EVOLVING IN A LOSSY 

GAUSSIAN CHANNEL 

In this section we test the effectiveness of our criteria, by 
applying them to typical quantum states that are of relevance 
to the quantum optics community. We shall consider pure, 
non-Gaussian states evolving in a lossy channel, and test their 
quantum non-Gaussianity after such evolution. Specifically, 



we focus on the family of quantum channels associated to the 
Markovian master equation 

dg . 7 



(a'ag - 



ga*a) 



(28) 



The resulting time evolution, characterized by the parameter 
models both the incoherent loss of photons 



1 



^7t 



in a dissipative zero temperature environment, and inefficient 
detectors with an efficiency parameter 77 = 1 — e. The evolved 
state £t{go) can be equivalently derived by considering the 
action of a beam splitter with reflectivity e, which couples the 
system to an ancillary mode prepared in a vacuum state. The 
corresponding average photon number reads 

n^ = Tr[£^{go)a^ a] = (1 - e) 



"•0 



(29) 



where uq — Tr[goa'a] is the initial average photon number. 
It is well known that, for e > 0.5 (i.e. for detector efficiencies 
7] < 0.5), no negativity of the Wigner function can be observed. 
We will focus then on the violation of our criteria for larger 
values of e, which ensures that the evolved states have a positive 
Wigner function. 

Notice that the quantum map £f^ is a Gaussian map. As a 
consequence, by combining the divisibility property of the 
map (inherited from the Markovian structure of Eq. (|28|) and 
Criterionl2] if a violation is observed for a given loss parameter 
e, then the state is quantum non-Gaussian for any lower value 
e < e 14 111 . For this reason we will focus on the maximum 
values of the loss parameter e for which a violation of the 
bounds is observed, i.e. 



eWje]=max{e : A,[£,{g)] < 0} , 

e^^l[g] ^ max{e : 3£g s.t. A2[£e{g),£G] < 0} 



(30) 
(31) 



In what follows, we start by focusing on Criterion [T] and thus 
we will look for negative values of the non-Gaussianity indi- 
cator Ai[g] defined in Eq. (26 1. We will consider different 
families of states, namely Fock states, photon-added coher- 
ent states and photon-subtracted squeezed states. In section 



IV B we will study how to improve the results obtained, by 
considering the second criterion and thus by studying the non- 
Gaussianity indicator A2 [g, £g] ■ 



A. Violation of the first criterion 

1. Fock states 

Let us start by considering Fock states |m), that is the eigen- 
states of the number operator: d''a\m) = m\m). A fock state 
evolved in a lossy channel can be written as a mixture of Fock 

states as 



Al[p] 
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FIG. 2. (Left) Non-Gaussianity indicator Ai[£'e(|m)(r7i|)] for the 
first three Fock states: red-dotted line: m — 1; green-dashed line: 
m = 2; blue-solid line: m = 3. 
(Right) Maximum value of the noise parameter e|„ax such that the 



bound 1 14 1 is violated for the state £^{\m){m\), as a function of the 
Fock number m. 



We recall here that the Wigner function at the origin is 
proportional to the expectation value of the parity operator 

n= (-)"*'', that is 

W[gm - ^Tv[gn] = ^ (Peven - ^odd) , (34) 

where Pgven (^odd) represents the probability of detecting an 



even (odd) number of photons. By using Eq. (32 1 one obtains 

W[£,{\m){m\m^-i2e~ir, (35) 

TT 

and thus the non-Gaussianity indicator reads 

(36) 

The behavior of Ai[E{\m){m\] as a function of e for the 
first three Fock states is plotted in Fig. |2](left panel). One 
can observe that the criterion works really well for the Fock 
state |1), which is proven to be quantum non-Gaussian for 
all values of e < 1. For the Fock states |2) and |3), a non 
monotonous behavior of Ai is observed as a function of the 
loss parameter Still, negative values of the non-Gaussian indi- 
cator are observed in the region of interest e > 0.5. However, 

the maximum value of the noise parameter efrax decreases 
monotonically as a function of m, as shown in Fig. l2] (right 
panel). By increasing the Fock number m, it settles to the 
asymptotic value efrax — ^ 0.5. As one would expect by looking 



at the bound in Eq. ( 14 1, for high values of the average photon 



number, the criterion becomes practically equivalent to the 
detection of negativity of the Wigner function, and thus the 
maximum noise corresponds to e = 0.5. 



with 



£,{\m}{m\) = J2^i,m{^)\l){l\ 



1=0 



ai,m{^)=[ ^ )(1-6)V 



(32) 



(33) 



2. Photon-added coherent states 

A photon-added coherent (PAC) state is defined as 

1 



IV'pac) 



yiTR 



aT \a) . 



(37) 
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FIG. 3. (Left) Non-Gaussianity indicator Ai[£:,(|i/'pac)(V'pacl)l for FIG. 4. (Left) Non-Gaussianity indicator Ai[£:,(j^pss){i/)pss|)] for 



PAC states as a function of e and for different values of a: red-dotted 
line: a — 0.2; green-dashed line: a = 0.4; blue-solid line: a = 0.6. 



bound 1 14 1 is violated for the state fe(|i/'pac){^pacl), as a function of 



the parameter a. 



PSS states as a function of e and for different values of r: red-dotted 
line: r — 0.1; green-dashed line: r — 0.3; blue-solid line: r = 0.5. 
(Right) Maximum value of the noise parameter e|„ax such that the 
bound l |14[ is violated for the state S^dippss) {ip pss\), as a function of 
the initial squeezing parameter r. 



The operation of photon-addition has been implemented in 
different contexts Il42l - l45l . and in particular non-Gaussianity 
and non-classicality of PAC states have been investigated in 



Being the non-Gaussianity indicator Ai [g] phase insensitive, 
we can consider a e R without loss of generality. The average 
photon number can be easily calculated obtaining 

/ / I t I / \ a^ + 3a^ 



;(pac) 



1 



1 



a^ 



while its Wigner function reads 

2 g-2(a-A)(Q-A-) 



W[\%,c)KX) 



-1- 



1 + a^ 
^2 + 4|AP 



2a(A + A*)) 



(38) 

(39) 
(40) 



The Wigner function of the state after the loss channel 
'£'e(|V'pac)(V'pac|) Can be evaluated by means of the formula 



W[£,ig)]i\)= d^X'K,{\,X')W[Q]i\'), (41) 



where 



2 2\\-\'^/l^^ 

K,{X,\') = —expl ^ 

7re e 



(42) 



The non-Gaussianity indicator Ai[f(:(|'!/)pac)(f/'pac|) can then 
be straightforwardly evaluated and is plotted in Fig. [3] (left) as 
a function of e for different values of a. We note that negative 
values of the indicator can be observed in an interval for the 
noise parameter e, which decreases with the increase of a. We 
can explain this feature by noting that, as a decreases, the 
PAC state approaches the Fock state |1): as a consequence 
its quantum non-Gaussianity can be more easily detected via 
Criterion [T| in particular due to the minimum value of the 
Wigner function approaching the origin of the phase space. We 
plotted in Fig. p|(right) the maximum value emax at which the 
violation of the bound is observed as a function of a. Similarly 
to Fock states, we observe that by increasing the energy this 
value tends to the asymptotic value emax —> 0.5. 



3. Photon-subtracted squeezed states 

Let us consider now another important class of non-Gaussian 
states that can be engineered with current technology. The 
photon-subtracted squeezed (PSS) states are defined as 



l^pss) = 



1 



sinh 



-aSir) |0) 



(43) 



For low values of squeezing, these states approximate the 
Schrodinger kitten states, that is, superpositions of coherent 
states I ± a) with opposite phase and small amplitude (|a| < 1) 
P6l . The generation of this kind of states has been demon- 
strated experimentally f4T-30l, and it relies on performing 
conditional photon number measurements. 
Without loss of generality we shall consider a real squeezing 
parameter r g R; the corresponding average photon number 
of a PSS state reads 



,(pss) 



3 sinh 



1 



(44) 



while its Wigner function is 



W^[|^pss)](A) = 



_ _^ -2|A|'^ cosh 2r+(A^+A*^) sinh 2r, 



[1 - 4| Ap cosh 2r + 2{\^ + A*^) sinh 2r] 



(45) 



As for the PAC states, the Wigner function of the evolved state 
can be evaluated by means of Eq. (J4Tll and the non-Gaussianity 
indicator Ai [£^{\il]ps^ (''/'pssl] can be evaluated accordingly. Its 
behavior as a function of e and for different values of the 
squeezing factor r is plotted in Fig. l4](left). In the right panel of 
Fig. |4|we plot the maximum noise parameter emax as a function 
of the squeezing parameter r, observing the same behavior 

obtained for Fock and PAC states: the value of emax decreases 
monotonically with the energy of the state, approaching the 

asymptotic value Cmax — > 0.5. 
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FIG. 5. Contour plot of the Wigner function of the photon-added 
coherent state j^pss) for a = 1. The minimum of the Wigner function 
is not at the origin of the phase-space, and the state has non-zero first 
moments. 



B. Violation of the second criterion 

We will now show how the second criterion, which is based 
on the violation of the inequality ( [ST) , can be exploited in order 
to improve the results shown in the previous section. Since in 
this case one can optimize the procedure over an additional 

tmax 



> eLx. The sim- 



Gaussian channel, in general one has €„ 

plest Gaussian maps that one can consider are displacement 

and squeezing operations; coiTespondingly we are going to 



seek violation of the bounds described by Eqs. ( |23| ) and (|24 1 
As anticipated in Sec 



III 



these new criteria are useful for 
states which are not phase invariant: the paradigmatic exam- 
ples are states displaced in the phase-space, that is, having 
the minimum of the Wigner function outside the origin, or 
states that exhibit squeezing in a certain quadrature. Due to 
this fact, the bounds based on Eqs. ( p3] l and ( |24| l cannot help 
in optimizing the results we obtained for Fock states. We will 
focus then on the other classes of states we introduced, that is 
PAC and PSS states. 



1. Photon-added coherent states 

By looking at the PAC state Wigner function in Fig. l5| one 
observes that its minimum is not at the origin of the phase 
space. Moreover, these states have non-zero first moments, 
implying that one can decrease their average photon number 
by applying an appropriate displacement. Both observations 
suggest that it is possible to decrease the value the quantum 
non-Gaussianity indicator defined in Eq. ([27]i, 



Apac(/3) = A2[f.(|V'pac)(^pac|),2?^] 



(46) 



by means of a displacement operation 2?g(g) = D{/3)gD{f3y . 
To evaluate Apac(/3) according to Eq. (46 1 one has simply to 




FIG. 6. (Left) Non-Gaussianity indicator Apac(/3) as a function of the 
additional displacement parameter f}, for e — 0.8 and for different 
values of the initial parameter a: : red-dotted line: a — 0.2; green- 
dashed line: a = 0.4; blue-solid line: a — 0.6. 
(Right) Optimized non-Gaussianity indicator Apac(/3opt) as a function 
of e and for different values of a, where the displacement parameter 
/3opt has been chosen as in Eq. l |49[ >: a — 0.2; green-dashed line: 
a = 0.4; blue-solid line: a = 0.6. 



evaluate the Wigner function of the state g — £'e(|V'pac) (V'pacl) 
in a displaced point in the phase space, i.e. W[g]{—f3), and its 
average photon number 

+ VT^X/3*{a)o + l3{a^)a) , (47) 

where {A)o = {^|Jo\A\^po) , and for I^Ao) = iV'pac), 



— I n < \ r^ — 



(a)o = (aVo = 



(48) 



Our goal is then to minimize Apac(/3) over the possible dis- 
placement parameters (3. 

In Fig. l6](left) we plot Apac(/3) as a function of (3 for different 
values of the coherent state parameter a and for e = 0.8. We 
observe that, while for /? = the bound is not always violated, 
it is possible to find values such that Apac(/3) < and thus 
prove that the state is quantum non-Gaussian. Unfortunately 
the optimal value /3opt, which minimizes Apac(/3), can not be 
obtained analytically. However we observed that for large 
values of e and for a > 1.5 one can approximate it as 



Popt 



yVi 



^,^-7t/2 



The behavior of Apac(/3opt) as a function of e shown in Fi g. [6| 
(right), for different values of a and fixing /3opt as in Eq. ( |49[ ). 
If we compare this with Fig. [3] not only we observe an im- 
provement in our capacity to witness quantum non-Gaussianity 
for these states, but we also see that Apac(/3opt) remains nega- 
tive for all values of e. Indeed, numerical investigations seem 
to suggest that emax — 1 for all the possible values of a: we 
indeed conjecture that any initial PAC state remains quantum 
non-Gaussian during the lossy evolution induced by Eq. ( |28] l, 
and that this feature can be captured by our second criterion. 
However, as one can observe from Fig. l6] (right), the non- 
Gaussianity indicator approaches zero quite fast with both a 
and e, and thus it may be more challenging to detect its neg- 
ativity in an actual experiment for states with a high average 
photon number and for large losses. 
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FIG. 7. Contour plot of the Wigner function of the photon-subtracted 
squeezed state jV'pss) for r = 0.3. The minimum of the Wigner 
function is at the origin of the phase-space, and the state exhibits 
squeezing in one of the quadratures. 
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FIG. 8. (Left) Non-Gaussianity indicator Apss(s) as a function of the 
additional squeezing parameter s, for e = 0.7 and for different values 
of the initial parameter r: red-dotted line: r = 0.1; green-dashed line: 
r = 0.3; blue-solid line: r = 0.5. 

(Right) Optimized non-Gaussianity indicator Apss(sopt) as a function 
of e and for different values of r, where the squeezing parameter is 
given by: red-dotted line: r = 0.1; green-dashed line: r — 0.3; 
blue-solid line: r = 0.5. 



The optimized quantum non-Gaussianity indicator Aps5(sopt) 



2. Photon-subtracted squeezed states 

Like PAC states inherit a displacement in phase space from 
the initial coherent states, PSS states inherit squeezing, as we 
can observe by looking at the Wigner function in Fig. |7] This 
motivates us to make use of Corollary l2j and thus optimize the 
non-Gaussianity indicator in Eq. ( p7| as 

Apss(s) - ^2%{\%ss){%sslS,] , (50) 

that is by considering an additional squeezing operation 

Ss{q) = 5(5) giS''''(s) on the evolved state g — i?(|V'ps5)(V'pss|)- 
As pointed out in the proof of inequality ( |24[ l, the Wigner 
function at the origin is invariant under squeezing operations. 
Hence, the optimal value Sopt that minimizes Apa5(s) coincides 
with the value which minimizes the average photon number of 
S, = S{s)gSHs), 



i(p==)(.«) = 



W = (i-6) 



,(pss) 






(51) 



where /it — cosh t, vt ~ sinh t and for an initial PSS state 
(with a real squeezing parameter r). 



(a^>0 = (a"^^)o = ilJLrVr 



(52) 



The behavior of Apas(s) as a function of the additional squeez- 
ing s is plotted in Fig. [8] As we observed in the previous case, 
the optimised criterion works in cases where the bound \\A\ 
(corresponding to s = 0) was not violated. 
Moreover the optimal squeezing value can be evaluated analyt- 
ically, yielding 



5opt 



Mopt 



— arccosh(/iopt) 





FIG. 9. Maximum values of the noise parameter efnax and emax, 
obtained respectively by means of the optimized and not-optimized 
criteria, for the state fe(j'(/'pss)(^pss|)> as a function of the initial 



(2) , 



.(1) 



squeezing parameter r. Red-dotted line: Emax; blue-solid line: e 

is plotted in Fig. |8] (right), where we observe that negative 
values are obtained for large values of losses. However, while 
for PAC states we had evidence that the maximum value of 
1 for all the possible initial states, this is no 



losses is e 



(2) 



(2) 



longer true for PSS states. The behavior of efrax ^s a function 
of r is plotted in Fig. l9] together with the previously obtained 

Emax- We can notice the big improvement in our detection 

(53) capability, obtained by exploiting Corollary |2l however for 

(2) 
large values of r we still observe that Cmax decreases towards 

(54) the same limiting value e 



(2) 



0.5. Moreover, as it can be 



observed in Fig. l8](right), the indicator Apa5(sopt) approaches 
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zero by increasing r and e, and thus also in this case it can 
become challenging to witness quantum non-Gaussianity with 
our methods, in experiments with large values of the initial 
squeezing r and large losses. 



V. CONCLUSIONS 

We have presented a set of criteria to detect quantum non- 
Gaussian states, that is, states that can not be expressed as 
mixtures of Gaussian states. The first criterion is based on 
seeking the violation of a lower bound for the values that 
the Wigner function can take at the origin, depending only 
on the average photon number of the state. To verify the 
effectiveness of the criterion, we considered the evolution of 
non-Gaussian pure states in a lossy Gaussian channel, looking 
for the maximum value of the noise where such bound is 
violated. We observed that the criterion works well, detecting 
quantum non-Gaussianity in the non-trivial region of the noise 
parameters where no negativity of the Wigner function can be 
observed. 

We have also shown how the criterion can be generalized 
and improved, by optimising over additional Gaussian opera- 
tions applied to the states of interest. Notice that in a possible 
experimental implementation one does not need to perform 
such additional Gaussian operations, such as displacement or 
squeezing, in the actual experiment. Indeed, it suffices to use 
the data obtained on the state itself, and then apply suitable 



post-processing to evaluate the optimized non-Gaussianity in- 
dicator 

Our criterion, which expresses a sufficient condition for 
quantum non-Gaussianity, shares some similarities with Hud- 
son's theorem for pure Gaussian states, in the sense that it 
establishes a relationship between the concept of Gaussianity 
(combined with classical mixing), and the possible values that 
a Wigner function can take. The successful implementation 
of our criteria corresponds to the measurement of the Wigner 
function at the origin of the phase space which, in turn, corre- 
sponds to the (photon) parity of the state under investigation. 
This may be obtained with current technology by direct parity 
measurement 1.51,1 . or by reconstruction of the photon distri- 
bution either by tomographic reconstruction or by the on/off 
method (5T-'65]. When the criterion is satisfied, one can con- 
firm that the quantum state at their disposal has been generated 
by means of a highly non-linear process, even in the cases 
where, perhaps due to inefficient detectors or other types of 
noise, negativity of the Wigner function can not be detected. 
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